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The basic equations in the theory of ionospheric cross modulation are reviewed. The 
suggestion by Rumi that the variations in the electron density, caused by perturbations of 
the attachment coefficient, can contribute to the total cross modulation is considered. It is 
found that the cross modulation resulting from these variations is negligible compared to 
the cross modulation resulting from the variations in the electron collision frequency in the 
region above about 40 km. In the 30 km region, however, the two components are approx- 
imately equal. The fractional change in electron energy as predicted by the original theory 
of cross modulation, introduced by Bailey and Marty n in 1934, is compared with the same- 
quantity as predicted by the alternate theory of cross modulation proposed by Huxley in 
1953. 

Cross-modulation profiles are presented for these two theories, corresponding to various 
model ionospheres, and discussed in light of previously published cross-modulation observa- 
tions from College, Alaska. It is concluded that neither of the theories hold over the entire 
D region and that a new theory of ionospheric cross modulation is necessary. 

The requirements of such a theory, in order to be consistent with observations, are the 
following: for given ionospheric conditions the cross modulation should change sign in 
approximately the same height region as is predicted by the original theory of cross modula- 
tion but the absolute magnitude of the cross modulation, at least in the lower D region, 
should be equal to or greater than the value predicted by the alternate 1 theory of cross modu- 
lation (which, in turn, is greater than the value predicted by the original theory). A recent 
theoretical investigation by Altshuler, which indicates an energy dependence for the frac- 
tional change in electron energy similar to the predictions of the original theory, seems to 
satisfy the first requirement. The second requirement remains to be investigated. 



1. Introduction 

The phenomenon of ionospheric cross modulation 
is commonly associated with Tellegen's [1933] ob- 
servation that the reception of a Beromiinster, 
Switzerland radio program (on 650 kc/s) was, at 
times, marred by the presence of a Luxembourg- 
radio program (on 252 kc/s) in the background. 
After eliminating all possible local causes, Tellegen 
suggested that the interference effect was due to the 
interaction of the two radio waves in the ionosphere. 

This phenomenon has been known under various 
titles: "Luxembourg effect/' "radio wave inter- 
action/ ' and "ionospheric cross modulation.' ' 
Huxley and Ratcliffe [1949], in their review article 
on ionospheric cross modulation, suggested that the 
latter name be adopted and that the interacting 
radio waves be identified as the "wanted' ' and "dis- 
turbing" waves. 

Bailey and Martyn [1934] explained the cross- 
modulation effect in terms of the absorption of the 
radio waves in the ionosphere. The absorption of 
a radio wave as it travels in an ionized medium is 
dependent on the electron density and on the elec- 
tron collision frequency. The radio wave from a 
high powered disturbing transmitter (200 kW in 
the case of the Luxembourg transmitter) literally 
heats up the free electrons in the ionosphere, which 
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results in an increased electron velocity and electron 
collision frequency. The increase in the electron 
collision frequency alters the absorption of the 
wanted wave as it passes through the disturbed 
region. If the disturbing wave is amplitude modu- 
lated, the energy transferred to the electrons, and 
hence the absorption of the wanted wave, will vary 
in accordance with the degree of amplitude modula- 
tion on the disturbing wave. Thus, a fraction of 
the amplitude modulation present on the disturbing 
wave is transferred to the wanted wave giving rise 
to ionospheric cross modulation. 

The original theory of Bailey and Martyn [1934] 
was revised by Bailey [1937] and was later put into 
a form more familiar to the ionospheric physicist 
[Shaw 1951, and references therein]. Huxley [1953] 
proposed an alternate development of the theory of 
ionospheric cross modulation in an attempt to rec- 
oncile inconsistencies between the apparent behavior 
of electrons in the laboratory and in the ionosphere. 
His alternate theory was rejected [Huxley, 1955; 
Fejer, 1955] on the basis that it was inconsistent 
with the magneto-ionic theory. Also, the cross- 
modulation observations conducted in Norway 
[Landmark and Lied, 1961; Holt, Landmark, and 
Lied, 1961; Barrington and Thrane, 1962] and in 
Alaska [Rumi, 1961; Flock and Benson, 1961; Rumi, 
1962a, 1962b] seem to be in agreement with the 
original theory. 

Some of the recent observations in Alaska, how- 
ever, do not agree with the predictions of the original 
theory of cross modulation [Benson, 1963]. Pre- 
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liminary calculations indicated that the cross modu- 
lation, observed below about 50 km were in better 
agreement with the alternate theory proposed by 
Huxley [1953] than with the original theory of Bailey 
and Martyn [1934]. The above work [Benson, 1963] 
considered only the original theory of cross modula- 
tion in detail, and it is the purpose of this paper to 
critically analyze and compare the above two theories 
of ionospheric cross modulation. The validity of 
the suggestion given by Rumi [1962a], indicating 
that the contribution to the total cross modulation 
caused by variations in the electron density can be 
comparable to the contribution caused by varia- 
tions in the electron collision frequency, will also be 
examined. 



2. Basic Equations 

Consider the situation where the wanted wave is 
traveling vertically downward and the disturbing 
wave is traveling vertically upward. The absorption 
of the wanted wave as it passes through an in- 
finitesimal homogeneous layer of thickness dh is 
given by 

E=E e- K * dh (1) 



where, referring to the wanted wave, 

E= amplitude of the electric field of the emergent 

wave 
E = amplitude ol the electric field of the incident 

wave 
K w = absorption coefficient. 



The absorption coe K is given by 

i 



K=(5.31X10- 8 )^/|e 6/2 (a)| 



when the absorption is of the nondeviative type, 
i.e., the real component of the refractive index is 
unity, and the radio wave propagation is parallel 
to the direction of the earth's magnetic field [Sen and 
Wyller, 1960, combining their equations 39, 41, 
and 43]. In (2), the symbols have the following 
meaning : 

N= electron density (m~ 3 ) 

p m =the mean electron collision frequency asso- 
ciated with the most probable electron 
speed (sec -1 ) 






de 



a= for the ordinary wave and - - for the 

extraordinary wave (radians). 
(*) = angular radiofrequency (radians/sec), and 
s= angular gyrofrequency of the electrons due 

to the earth's total magnetic field in the 

region of interest (radians/sec). 



The script C integrals & p (a) have been tabulated, 
for a ranging from zero to 20, by Dingle, Arnt, and 
Rov [1957]. Rationalized inks units were used in 

(2): 

Assume that the downward traveling wanted 
wave passes through this infinitesimal layer after 
the passage of the upward traveling disturbing wave. 
The energy absorbed from the disturbing wave in- 
creases the electron collision frequency v m which 
causes the amplitude of the wanted wave to change 
by an amount 

aW^) (A,„,)„ 
\Ov m / h > 



-E 






(Av m ) h -dh' 



(3) 



where h f is the height of the disturbed layer of thick- 
ness dh r . 

Rumi [1962a] also considered the contribution to 
the cross modulation caused by variations in the 
electron density N. These variations in N are at- 
tributed to perturbations of the attachment proc- 
esses rather than to direct ionization. His analysis 
is based on recent laboratory measurements of the 
attachment of slow electrons in oxygen which indi- 
cate that the electron attachment is a three-body 
process with an attachment coefficient that increases 
with increasing electron energy in the energy range 
appropriate to the ionospheric I) region [Chanin, 
Phelps, and Biondi, 1959]. Thus, an increase in 
the electron energy, caused by the absorption of 
energy from the disturbing wave, should produce a 
decrease in the electron density. 

When the above effect is considered, (3) becomes 



AE-- 






i -+w 4 



}. 



dh'. 



(4) 



If the wanted and disturbing waves are both pulse 
modulated, and if the pulse repetition rate of the 
disturbing transmitter is one half that of the wanted 
transmitter, then only every other wanted pulse will 
be altered in amplitude (provided that the repetition 
period of the wanted pulses is chosen to be much 
longer than the deca}^ time of the excess electron 
temperature caused by the disturbing pulse). In 
this case the cross modulation is defined as the frac- 
tional change in echo amplitude, i.e., AE/E where E 
is the amplitude of one received echo pulse and 
E— AE is the amplitude of the adjacent pulse. Thus 
the cross modulation resulting from one infinitesimal 
disturbed layer at the height W is given by 



AT- 



AE 
' E 



dKu AN\ dh'. 

J h' 



(5) 



If the wanted wave is a continuous wave (such as a 
satellite signal or cosmic noise) the same comments 
apply provided that the wanted wave is sampled at 
a repetition rate that is twice that of the disturbing 
pulse-modulated transmitter. 
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Equation (5) gives the cross modulation resulting 
from only one infinitesimal disturbed layer. The 
downgoing wanted wave will encounter a succession 
of such disturbed layers below the height h where 
it first encounters the upgoing disturbing pulse. 
The total cross modulation observed on the received 
wanted wave is given bv 



T--=J^AT= 



- J^A Vm+ ^ANy dh>. 



(6) 



The evaluation of the partial derivatives that 
enter in (6) is straight forward when the longitudinal 
expression for K W) as given in (2), is used. This 
compact expression can be retained, rather than 
using the more involved expression for arbitrary di- 
rection of propagation, even when the direction of 
propagation deviates from the direction of the 
earth's magnetic field if s is replaced by u L =|s cos <j>\, 
where </> is the angle between the propagation vector 
and the earth's lines of magnetic force. Just as in 
the classical magneto-ionic theory, the above " quasi- 
longitudinal" approximation in the generalized 
magneto-ionic theory must satisfy certain conditions 
before it can be used [Benson, 1964a]. In the 1-4 
Mc/s frequency range, common to many experiments 
in ionospheric cross modulation, greater accuracy is 
obtained by using the unmodified longitudinal equa- 
tion (2), even when (/>^(), rather than by introducing 
o) L . In the cross modulation experiment in Alaska, 
the lowest frequency used is 4.865 Mc/s and 0=13°; 
under these conditions the quasi-longitudinal ex- 
pression for K w , i.e., using o) L in place of s in (2), can 
be used with an accuracy of better than 1 percent 
(with respect to the more general expression for 
arbitrary value of </>). The derivatives of K w , 
required in (6), are then given by 



dVm 



-(5.31X10- 6 ) 



where 



^=(5.31X10- 6 )^|(i 5/2 (a)| 



F(a)=|(S 5/2 («)-^6 7/2 (a). 



(7) 



(S) 



In deriving (7), the dummy variable e (the normal- 
ized electron energy) in the script C integral was 
written in its complete form as e = Q/kd, where Q is 
the energy of any given electron (which is the true 
dummy variable) and kd is the thermal energy 
associated with the electron having the most probable 
speed of all the electrons at the temperature d{°K). 
The assumption d=bv nl , where b is a constant, was 
then used in order to carry out the differentiation. 
This assumption, that the electron collision frequency 
is directly proportional to the electron energy, was 
also used by Sen and Wyller [1960] to obtain their 
generalized magneto-ionic equations, and it is based 
on recent laboratory studies with N 2 by Phelps and 



Pack [1959]. Using the above assumption, the 
Av m term can be expressed as 



Av m =v m (Ad/6) 



(9) 



w here (A0/0) is the fractional increase in the electron 
temperature caused by the energy absorbed from the 
disturbing wave. The derivation of a similar ex- 
pression for the AN term, and an evaluation of the 
importance of the AN contribution, will be considered 
next. 

3. Variations in the Electron Density 

In order to obtain an expression relating the term 
AN to Ad the following equations, which describe 
the rate of change of charge density in the ionosphere, 
will be considered: 



dt 



: q.+7pN~ +y c nN~ —a,,NN+ 



- 0[ n(0 2 ) )N- k[ n(0 2 ) ]W (1 0) 



ftN ~ = -y N-- 
dt 7 " 


-y> l N--a i A T -N+ 




+p[n(0 2 )]N+K[n(0 2 )] 2 N (11) 


N+=N+N~ 


(12) 



where 

N= electron density 
AT" = negative ion density 
iV + =positive ion density 
[n (02)]= number density of molecular oxygen 
n= number density of neutral particles 
g= electron production rate 
a d = electron-ion dissociative recombination co- 
efficient (effective value) 
a, = ion-ion recombination coefficient 
|S = 2 bod} 7 attachment coefficient 
k=3 body attachment coefficient 
7 C = collisional detachment coefficient 
y p = photodetachment coefficient. 

Rumi [1962a] considered two cases of the equilib- 
rium equation (which follows from (10), (11), and 

, 10 . , dN dN~ A 
(12) when —-==-^-=0 J 



namelv, 



dt dt 

( L =(l+X)(a d + \a i )N 2 



(13) 



and 



N=l (£) when X»l and \a.»a d (14) 
X \oti/ 

N=(-^t) when x »l and Xa / «a d (15) 
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where \=N~/N. From (14), which applies to the 
lower D region, he derives the following expression 

for AiV: 



AN= 



and presents a curve of 6 — — — versus which is 

based on the data of Chanin, Phelps, and Biondi 
[1959]. 

In the above analysis, however, it is inherently 
assumed that the time constant for the variations in 
the electron density is the same as the time constant 
for the variations in the electron collision frequency. 
A recent report by Molmud, Altshuler, and Gardner 
[1962] indicates that this is not the case in the iono- 
spheric D region. The above authors discuss the 
time constant for electron density changes resulting 
from electron energy changes in their study of a 
method to reduce the electron density in the iono- 
spheric D region by means of high-powered ground- 
based transmitters. Their study is based on the 
fact that the time constant for electron density 
variations is long compared to the time constant for 
average electron energy variations. 

Molmud, Altshuler, and Gardner [1962] combined 
(10), (11), and (12) by eliminating N~ while retaining 
the time parameter t. They then consider the solu- 
tion of the resulting differential equation under the 
following initial conditions : for t<C0, = 6 and steady 
state conditions prevail; for t>0, 0=0i. The rise in 
electron temperature, however, is generally not 
instantaneous, even if a pulse modulated disturbing 
wave is used. (This subject is discussed in detail in 
the next section.) Thus the solution of the differ- 
ential equation for N, under the above initial con- 
ditions, will be useful only in determing an upper 
limit for the variations in N due to the true variations 
in 6 caused by the disturbing pulse. Combining 
(10), (11), and (12), using the method of Molmud, 
Altshuler, and Gardner [1962], gives 



2^=A-2BN 
dt 



(16) 



where 



A=a + 2{y v +y c n)N++a i {N + ) 2 
B=y p +y e n+(^^ N++k[u(0 2 )Y. 

In deriving (16) it has been assumed that the changes 
in N + are slow compared to the changes in N. 
Molmud, Altshuler, and Gardner [1962] show that 
this is a reasonable assumption. 

The rate coefficients, that appear in (16), will be 
assigned the following values [Crain 1961]: 

a rf =1.0X10" 7 cm 3 /sec 

a t =6[p (atm)]r 5/2 +10- 9 cm 3 /sec 



7 P =0 during the night 

=0.4 sec -1 during the day 
Yc^<100 p (atm) sec -1 

q m =(a d +hx i )NN + (the value for equilibrium con- 
ditions). 

The positive ion density is given by N + =(l-{-\)N, 
and X=iV"7iVcan be approximated by 

^=[n{0 2 )] 2 K/(y p +y c n) [Rumi, eq (13), 1962a]. 

The following approximate expression for the three 
body attachment coefficient k can be derived by ex- 
trapolating the 2 curve in figure 3 of the paper by 
Chanin, Phelps, and Biondi [1959]: 

k (cm 6 /sec) ^ -2.20 X 10- 30 + (1.63 X 10- 32 )<9. 

Sample calculations using the above expressions 
indicate that the parameter A, in (16), is dominated 
by the term 2(7 7 ,+Y c ^)iV + , which is independent of 
the electron energy. Thus, to first approximation, 
A(e )=A(di)==A. The parameter B, on the other 
hand, is dominated by the energy dependent term 
k[u(0 2 )] 2 . The initial conditions appropriate to (16) 
are then given by 



B=B 

N=N j 



H<0 




B=B X 

N=N. 



► *>o 



where 2Vo is the equilibrium value of iVin the absence 
of the disturbing pulse (2<0) which is found from 
(16), with dN/dt=Q, to be 

N =A/2B . (17) 

The integration of (16) then proceeds as follows: 

2dN 



J N o 



A-2B X N 



f dt 



K™>=-*" 



Combining (17) and (18) gives 

=2§ 2 (Bo- Br) (l-e- Bt ) 



(18) 



(19) 



where the distinction between B and B x is made 
only in the difference term (Bq—B^), since the differ- 
ence between the two quantities is slight (the maxi- 
mum value of Ad/9 is of the order of a few percent 
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or less). Consider the term (B —Bi). From (16), 
and the values for the rate coefficients given below 
(16), (Bq—Bx) can be written as 

B Q ^B 1 =lJ2[(a i ) -(a i ) 1 ]m+(K Q --K 1 )[n(0 2 )] 2 

tfl5/2_ ^5/2-1 

+ (1.63XlO- 32 )(0 o -0 1 )[w(O 2 )] !! . 

Let6>i=0 o +A0, then, 

/ A0V /2 / Ad\ 



9 e ) 



and 



Thus 



j}/»_jg/»~5/2 (y) 5/2 . 



Hi 



5[p(atm)]iV 4 



^ |i , v ,^y J .. _ (L63X10 - 3 2)^(0 2 )] 2 



Ad 



and (19) becomes 



AN~H(h') ^ 

C7 



(20) 



where 



fl r (A / ) = 



A (15 [p(atm)]iV+ 
~2B 2 \ 2 5/2 



-(1.63X10- 32 )<9M0 2 )] 2 j (1— e- BI ) cm" 3 

and A and B are given in (16). Note: cgs units were 
used exclusively in the above equations because the 
rate coefficients are most commonly expressed in 
these units. In the following discussion, H(m~ 3 ) 
= 10 6 H (cm -3 ) will be used. 

Combining (6), (7), (8), (9), and (20) gives the 
following expression for the total cross modulation 
observed on the received wanted wave : 

T=^(5.31X10- 6 ) r {N(m-*)Y(a) 
^ Jo 

-H(m-^ i/2 (a)}-^dh'. (21) 

In (21), the term NY (a) determines the component 
of cross modulation due to variations in v m {CMv), 
and the term HG£y 2 (a) determines the component 
of cross modulation due to variations in N (CMN). 
From this equation, the importance of CMN, as 
compared with CMv, can be estimated. Sample 
calculations, using the model atmosphere given by 
Miller [1957] together with a model profile for the 
electron density [see Benson, 1963, profile N^], in- 
dicate that CMN is negligible compared with CMv 
except in the region below about 40 km. The results 
of these calculations are listed in the following table. 









CMN/CMv* 




h (km) 














u- 


= 17.5 


Mc/s 


/»= 


= 4.865 Mc/s 


30 




1. 






0. 9 


40 




. 1 






. 07 


50 




. 07 






. 009 



*CMN=—H(M-*)£ bn (a); CMv = N(m~*) F(o). 

The above figures are maximum figures appropriate 
to a disturbing pulse of 50 /xsec duration; i.e., in 
(20), t was set equal to 50X10 -6 sec. It must be 
kept in mind that the true contribution of CMN is 
less than is indicated by the above figures, especially 
in the region above about 40 km, because the in- 
crease in electron temperature is not instantaneous. 
(In the region below about 40 km, the increase is 
nearly instantaneous when the original theory of 
cross modulation is used and the above figures ap- 
proximate the true situation; the increase is not 
instantaneous, however, when the alternate theory 
of cross modulation is used, and in this case the true 
figures are less than is indicated above.) The above 
results are in agreement with the results of Rumi 
[1962a] in the lower ionosphere (30 km region) when 
the original theory of cross modulation is used, 
namely, that CMN is comparable to CMv, but 
differ greatly from his results above this region. 
This situation exists because the time constant 
{I IB), associated with AN in (20), increases with 
increasing altitude, which prevents significant 
changes in N from occurring during the duration 
of the disturbing pulse. Rumi's treatment did not 
include the time constant effect since the time 
parameter was eliminated from (10) and (11). 



4. The Term Ad/d 

The terms A^ and AiV, that enter in the cross 
modulation equation (6), have each been expressed 
in terms of the fractional increase in the electron 
temperature Ad/d due to the disturbing pulse [see (9) 
and (20)]. This term, Ad/6, is equivalent to AQ/Q, 
where Q, the thermal energy of the electron, is given 
by Q=akd (k is Boltzmann's constant and a is a 
constant which takes the value 1, 4/tt, or 3/2 depend- 
ing on whether the electron thermal energy is con- 
sidered to be associated with the most probable 
electron velocity v m , the mean electron velocity v, or 
the root-mean -square electron velocity v rma ). The 
value a= 1 will be used throughout this paper so as 
to be consistent with the use of v m in (2). The 
variation of Q is determined by the equation 



dQ. 
dt 



=w—v m 8Q 



(22) 



where w is the power supplied to each electron and 
5Q is the mean energy lost by an electron per collision 
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with a heavy particle. In the original theory of 
ionospheric cross modulation, Bailey and Martyn 
[1934] assumed that 



8Q=G(Q-Q ) 



(23) 



where G is a dimensionless constant generally 
referred to as the energy loss coefficient, Q is the 
thermal energy of an average electron after the 
passage of the disturbing wave, and Q is its thermal 
energy in the absence of the disturbing wave. 
Equation (22) then becomes 



dQ 
dt = 



--w—Gv m (Q—Q ). 



Since Q is considered to be constant, this equation 
can be written as 



d(AQ)_ 
dt 



~-w-Gv m {AQ) 



(24) 



where AQ=Q— Q . 

Consider the solution of (24) appropriate to a 
pulse modulated disturbing wave traveling vertically 
upward with a velocity c (assumed to be equal to 
the velocity of light in vacuum). Since the wanted 
wave is traveling vertically downward at the same 
velocity, any portion of the wanted wave is only in 
contact with the disturbing pulse for a time r/2 
where r is the duration of the disturbing pulse in 
seconds. Thus, if a portion of the wanted wave 
first encounters the leading edge of the disturbing 
pulse at the height h , it will encounter the trailing 
edge of the disturbing pulse at the height h —cr/2; 
accordingly, its position at any later time t is given 
by h —ct/2. 

The value of AQ at a height h' when h >h'>h — 
ct/2 (i.e., inside the pulse) can be obtained by 
integrating (24). In this integration, v m will be 
approximated by an average value v a appropriate 
to the region between h and h —cr/2. The integra- 
tion proceeds as follows: 



/. 



AQ={AQ) h , 



d(AQ) 



w-6v a (.AQ) 



T 

pP-gy.(A®v T 



!(»„-*'> /« 



dt 



2%-h' 






w 



-2(ho-h') 



CT 



Gv a ~ 



(25) 



when h Q >h'>h — 

The value of AQ at a height h' when /*/</k — cr/2, 
i.e., outside the pulse, can be obtained by integrating 
(24) with w=0. The equation then becomes 



x 



<*»*' d(AQ)_ 



■Qvm£' 



= 2(ft n -c T /2-&') lc 



dt' 



where (AQ) max is the value of AQ at A / =A — ct/2 



and t' is the time measured with reference to the 
passage of the trailing edge of the disturbing pulse. 
Since the decay of AQ at a given height is under 
consideration, v m is a constant and has been taken 
outside of the integral sign. Integrating the above 
equation gives 



(AQ) v = (AQ) max 6 



-2(/i -Cr/2-/O 



(h,„ 



(26) 



when /i/<^ — cr/2. 



Huxley [1953] proposed an alternate expression 
for the mean energy lost by an electron per collision 
when laboratory studies did not seem to confirm 
the assumption in (23). The development of his 
alternate expression was also encouraged by the 
discrepancy existing between the laboratory measure- 
ments of G and the ionospheric cross modulation 
measurements of the quantity Gv. Following 
Huxley [1953], the following expression can be 
derived : 

SQ=B ( -^ (27) 

where 5=3.18X10~ 23 j. The details leading to this 
expression are given in the appendix. An important 
point to notice is that (27) is based on a steady 
electric field, and thus only applies to the region 
where v m ^>^>j dy where j d is the frequency of the 
disturbing transmitter. Substituting (27) into (22) 
and proceeding in the same manner as was used to 
derive (25) and (26) gives 



^■^y^ r^y^'' ] 



(28) 



when h >h' >h —- 
and 

(A<2)»< = 



^Vmax 



r 2Bv m (h -CT/2-h') -\ 

"A^max L C Q2 J 



(29) 



when h'<Ji — ct/2 

where AQ max is the value of AQ from (28) when 
h'=h — ct/2. 

Sample calculations were performed in order to 
compare the values of AQ/Q obtained from the 
original theory (25) and (26) with the values of 
AQIQ obtained from the alternate theory (28) and 
(29). The results are presented in figure 1. In 
performing these calculations, the following value 
was used for w: 

P d gK d - P' aiz*F 
V77/7V2 e Jo 



w= 



where 



2irN{h') 2 



watts 



(30) 



P<z=peak pulse power radiated by the disturbing 
transmitter (watts) 
g= disturbing antenna gain factor 
K d — absorption coefficient for the disturbing wave 
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2.2X10 




20 

(h -h') km 



Figure 1. The fractional increase in the electron energy Q, caused by a rectangular shaped upward traveling disturbing pulse of 
50 ii sec duration, as predicted by the original theory of cross modulation [Bailey and Martyn, 1934] and the alternate theory 
of cross modulation proposed by Huxley [1953]. 

Each curve corresponds to a particular value of ho, the height where the leading edge of the upward traveling disturbing pulse first encounters the downwind I raveling 
wanted wave. The curves correspond to the following parameters: Pj=l()kW, ^5,/d=4.865 Mc/s (E),/ w = 17.5 Mc/s. 



(see appendix). Model electron density and collision 
frequency profiles of a previous publication were 
used in the calculations (A^ and v l of Benson [1963]). 
The parameters were chosen so as to agree with the 
recent cross-modulation experiment in Alaska which 
uses cosmic noise at a frequency of 17.5 Mc/s as the 
wanted wave [Benson, 1962]; i.e., 1\=\0 kW, 
#~5, and/ (/ = 4.865 Mc/s. The calculations corre- 
spond to the extraordinary component of the dis- 
turbing wave. Each curve in figure 1 indicates the 
fractional variations in electron energy AQJQ that 
are encountered by a portion of the downward 
traveling wanted wave as it passes the upward 
traveling disturbing pulse when the leading edge 
of the disturbing pulse is at a given height A - The 
curves are drawn for the case r=50 /zsec; i.e., the 
trailing edge of the disturbing pulse is located at the 
height h — crj2 = h — 7.5 km. 

Figure 1 indicates that the expected value for 
AQIQ is larger and that it decays much slower in 
the region outside (below) the disturbing pulse, for 
any given value of h , when the alternate theory is 
used rather than the original theory. The maximum 
value for AQjQ increases with increasing h up until 
h Q =50 km; the maximum value for A >50 km 
decreases with increasing h due to the absorption 
of energy from the disturbing wave in the lower 
regions and to the slower rise of AQ/Q in the upper 



regions. The discontinuity in (he curves labeled 
/l =60 km and h =7() km at the height h/ = h {) —7.5 
km, in both theories, results from (he approximation 
v m =v a inside the pulse (v a was evaluated at the 
height h'=h — 3 km). (The larger values of AQJQ 
outside the pulse are due to larger values of (AQJQ) umx 
in the lower regions.) The curves for A =60 km 
and h = 70 km do not satisfy the condition v m ^>^>f d , 
which must be satisfied in the alternate theory; this 
point will be considered further in the next section 
where the original and alternate theories are com- 
pared b}^ calculating the cross modulation expected 
for each theory for various model ionospheres. 

5. Model Cross-Modulation Profiles 

When the term AA 7 is neglected and the term 
Av m is written as v m (A6/6) from (9), then the equation 
for the total cross modulation (6) becomes 



--f(£>->0*- 



(31) 



The partial derivative in this equation is given by 
(7); the Ad/6 term is given by (25) and (26) when the 
original theory of cross modulation is used, and by 
(28) and (29) when Huxley's alternate theory of 
cross modulation is used. 
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The cross modulation was calculated from (31), 
with the aid of the IBM 1620 computer at the Uni- 
versity of Alaska, for several model ionospheres. 
In these calculations the script C integrals were 
evaluated using the tabulated values, for the argu- 
ment ranging from to 20, as given by Dingle, Arndt, 
and Roy [1957]. These integrals were approximated 
as 1/a 2 for values of the argument a greater than 20. 
The electron temperature 6 was assumed to be com- 
parable to the gas temperature, and the Fort Churchill, 
Canada, winter average rocket measurements of the 
D region gas temperature by Stroud, Nordberg, 
Bandeen, Bartman, and Titus [I960, see fig. 13] were 
used in the calculations. The electron collision 
frequency v m was evaluated from the equation v m = 
8.40 X10 7 p (mmHg) using the atmospheric pressure 
values for November 1956 above Fort Churchill as 
given by Lagow, Horowitz, and Ainsworth [1960, 
see fig. 5]. This equation is based on the laboratory 
studies in N 2 by Pack and Phelps [1961] and is dis- 
cussed in the following paper [Benson, 1964b]. 

Several model profiles for the electron density N 
were used in the calculations; these profiles are 
shown in figure 2. The curve labeled iVi (the curve 
without the bumps) represents a possible electron 
density profile in the arctic regions during relatively 
quiet ionospheric conditions. Each of the curves 
labeled N 2 , N z , N±, and iV 5 represent particular 
perturbations of the main curve iVi. These per- 
turbed curves should be thought of as short period 
electron density profiles resulting from transient 
irregularities in the D region electron density rather 
than as stable electron density profiles. The non- 
deviative absorption of a radio wave at a frequency 
of 27.6 Mc/s (a standard riometer frequency at 
College, Alaska) for each of the above electron 
density profiles is given in table 1. 




Figuee 2. These curves are not the result of an experimental 
program; they are merely model electron density profiles used 
in calculating the model cross-modulation profiles. 

The curves N2, Nz, JV4, and iVj, which represent particular perturbations of 
the main curve N\, are introduced for the sake of discussing the equations 
of cross modulation. 



Table 1. Radio wave absorption figures for the model electron 
density profiles 



Electron density profile 


Ni 


N 2 


Nz 


Ni 


N 5 


Absorption in decibels at 27.6 Mc/s 


1.1 


5.2 


1.2 


2.3 


1.9 







In the original theory of cross modulation, it is 
necessary to assign a value to the energy loss coeffi- 
cient G, which appears in (25) and (26), before 
model cross-modulation profiles can be computed. 
The value 6r=lX10~ 3 , which is consistent with 
laboratory measurements in air [Crompton, Huxley, 
and Sutton, 1953], was used in the present calcula- 
tions. 

In the calculations, the parameters of frequency 
and power were chosen so as to be consistent with 
the cross-modulation experiment at College, Alaska. 
Three separate conditions were considered: the first 
condition corresponds to the recent cross-modulation 
experiment which used cosmic noise as the wanted 
wave [Benson, 1962], the second condition corre- 
sponds to the original cross-modulation experiment 
which used a reflected pulse modulated radio wave 
for the wanted wave [Rumi, 1961]; and the third 
condition corresponds to the cross-modulation ex- 
periment that will be under way when the present 
construction of a new cosmic noise receiving antenna 
is completed. These parameters are presented in 
table 2. 

Table 2. Frequency and power parameters employed in the 
calculation of model cross-modulation profiles 





Disturbing transmitter 


Frequency 




Fre- 
quency a 


Power 


Antenna 
gain 
factor 


of wanted 
wave a 


Condition 1 

Condition 2 

Condition 3 


Mc/s 

4. 865(E) 
17.5 
17.5 


kW 

b20 

100 

100 


5 
150 
150 


Mc/s 
17.5 

4. 865(0) 
18.0 



a (E) designates extraordinary component; (O) designates ordinary component, 
b In the early phases of the experiment a value of 10 kW was used; the curves 
in figures 3 and 4 correspond to the value of 20 kW. 

The radiofrequencies listed in the above table are 
considerably higher than the frequencies used in other 
cross-modulation experiments (usually in the 1 to 3 
Mc/s range). One advantage of the higher frequen- 
cies is that the cross-modulation experiment is most 
efficient during periods of ionospheric disturbance 
which are common at College, Alaska; one disad- 
vantage of the higher frequencies is that the sensi- 
tivity of the experiment is decreased in the upper D 
region. The latter point is illustrated in figure 3 
where a large electron density perturbation in the 
upper D region (see fig. 2, profile iV 2 ) is seen to pro- 
duce only an insignificant perturbation in the cross- 
modulation profile. The cross-modulation profiles 
presented in figure 3 are based on a wanted frequency 
of 17.5 Mc/s (condition 1 of table 2). The sensitivity 
of the experiment to electron density perturbations 
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40 50 60 70 

HEIGHT, km 



Figure 3. A comparison of the expected cross modulation for 

Ni and N 2 of figure 2 when the original theory of cross modula- 
tion is employed, using the parameters of condition 1. 



in the upper D region increases considerably for lower 
frequencies. It is important to notice, however, that 
the lack of sensitivity indicated in figure 3 is partly 
due to the slow rise of AQ/Q in the upper D region 
(see fig. 1). 

In figures 4, 5, and 6, the ionosphierc cross- 
modulation profiles, corresponding to the model 
electron density profiles Ni, iV 3 , iV 4 , and N 5 , are 
presented for conditions 1, 2, and 3 (see table 2), 
respectively. Note: Because the absolute value of 
win (30) was used, the curves indicate a sign reversal 
of T was respect to the curves presented in an earlier 
publication [Benson, 1963]. In each figure the cross 
modulation predicted by the original theory is 
compared with the cross modulation predicted by 
the alternate theory up to an elevation of 55 km. 
The dashed curves, corresponding to the alternate 
theory, are terminated at 55 km because above this 
height (approximately) the radiofrequency of the 
disturbing transmitter fails to satisfy the condition 
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Figure 4. Cross-modulation profiles based on condition 1 of table 2. 
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Figure 5. Cross-modulation profiles based on condition 2 of table 2. 



fd<^<^v m , which is implied in the derivation of the 
alternate theory [see discussion following (27)]. Also, 
calculations based on the alternate theory predict 
that a significant portion of the excess electron 
energy AQ, caused by a given disturbing pulse, 
persists until the next sample of the wanted wave is 
taken when h is greater than about 55 km. 

At first glance it appears that the choice of frequency 
and power parameters as given in condition 2 is 
preferable to the choice of these parameters as given 
in conditions 1 and 3. All of the curves, however, 
correspond to the case when a usable wanted signal 
is available at the receiving antenna, and in condition 
2 the wanted signal, which is a radio wave at a 
frequency of 4.865 Mc/s that has been reflected from 
the E or F region of the ionosphere, is often com- 



pletely absorbed before reaching this antenna. For 
example, the total nondeviative absorption of the 
ordinary component of a reflected radio wave (two 
way path through the D region) at a frequency of 
4.865 Mc/s is of the order of 30 dB when the electron 
density profiles iVi, N 3 , N 4 , or N 5 are used. The 
absorption of the cosmic noise signal, which is the 
wanted wave in conditions 1 and 3, is much less than 
the absorption of the reflected radio wave in condi- 
tion 2 because the received cosmic noise is at a 
relatively high frequency and it travels through the 
D region only once. For example, the total non- 
deviative absorption for the cosmic noise at 17.5 
Mc/s (condition 1) and at 18.0 Mc/s (condition 3) is 
of the order of 3 dB or less when the electron density 
profiles Ni, N 3 , N 4j or N 5 are used. It is to be empha- 
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Figure 6. Cross-modulation profiles based on condition 3 of table 2. 
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sized that the above electron density profiles are 
merely model profiles used for the sake of comparing 
the two theories of cross modulation under discussion. 
During typical periods of ionospheric absorption at 
College, Alaska the absorption of radio waves is more 
severe than the absorption indicated by the above 
electron density profiles; under such conditions the 
echo technique (condition 2) becomes ineffective due 
to the excessive absorption of the wanted wave, and 
the cosmic noise cross-modulation technique is more 
advantageous [Benson, 1962]. 



The main points to observe in figures 4, 5, and 6 
are the following: 

(1) In all cases, T alt >T orig where T alt refers to 
the cross modulation predicted from the alternate 
theory and T ri g refers to the cross modulation pre- 
dicted from the original theory. 

(2) The ratio T &n /T 0Tlg is greatest when excessive 
ionization is present in the lower regions, i.e., for the 
electron density profiles iV 3 , N 4 , and N 5 . This is to 
be expected since the decay of AQ is slower, and thus 
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the cross-modulation contribution of the lower regions 
persists longer, when the alternate theory is used. 

(3) The ratio T^i t JT 0Tig is greater for condition 1 
than it is for conditions 2 and 3. This is to be ex- 
pected since the power term w for the disturbing 
wave, as given by (30), is different in condition 1 
than it is in conditions 2 and 3 and this term enters 
as the first power in the original theory [see (25)] 
but enters as the square root in both the coefficient 
and the argument of the hyperbolic tangent term of 
the alternate theory [see (28)]. 

(4) Significant cross modulation in the lower re- 
gions is predicted only by the alternate theory when 
the disturbing power is small (condition 1). 

(5) For a given value of the wanted frequency, 
the point where the cross modulation T changes 
from positive to negative depends only slightly on 
the form of the electron density profile when the 
original theory of cross modulation is used. For 
example, this "crossover" point is always between 
51 and 53 km for all of the model N profiles in con- 
ditions 1 and 3 where/^^17.5 Mc/s, and 18Mc/s re- 
spectively, and it is always between 60 and 61 km in 
condition 2, where f w =4.865 Mc/s (ordinary com- 
ponent). 

6. Discussion 

Figures 4, 5, and 6 indicate that the cross modula- 
tion predicted by the alternate theory [Huxley, 
1953] does not differ greatly from the cross modula- 
tion predicted by the original theory [Bailey and 
Martyn, 1934] when the D region electron density 
corresponds to the fairly moderate profile (for the 
arctic regions) given by Ni of figure 2. This is 
especially true when a disturbing transmitter of 
high power is used (see the cross-modulation curves 
corresponding to Ni in figs. 5 and 6). This is in 
marked contrast with the criticism of the alternate 
theory given by Huxley [1955] and Fejer [1955]. 
The criticism of Huxley [1955] was based on the cross- 
modulation measurements with obliquely traveling 
radio waves where the region of interaction was 
approximately 90 km. In these experiments fd~v m , 
whereas in the present discussion, concerned with the 
lower D region, f d <C<C.v m . Also, the classical 
magneto-ionic theory was used in the earlier ex- 
periments; the generalized magneto-ionic theory 
[Sen and Wyller, 1960] was used in the present 
calculations. 

The criticism of Fejer [1955] was based on cross- 
modulation observations which were confined to the 
upper D region where the condition f d <<^m, as 
required by the alternate theory, was not satisfied. 
Also, the classical magneto-ionic theory was used in 
the cross-modulation equation. 

From the above comments it appears that the 
objections to the alternate theory are not sufficient 
to completely neglect this theory in the lower D 
region. 

The observations of considerable cross modulation 
in the lower D region above College, Alaska are in 
better agreement with the alternate theory than with 
the original theory. For example, in the cosmic 
noise cross-modulation experiment at College, with 



the system parameters corresponding to condition 1, 
cross modulation was frequently observed in the 
40 to 45 km region [Benson, 1962 and 1963]; in the 
echo-type cross-modulation experiment at College, 
with the system parameters corresponding to con- 
dition 2, short duration periods of cross modulation 
were observed as low as 30 km [Rumi, 1961 and 
1962b; Flock and Benson, 1961]. The model elec- 
tron density profiles 2V 3 and N 4 were constructed to 
approximately simulate the former case and the 
profile N 5 was constructed to approximately simulate 
the latter case. The corresponding cross-modula- 
tion profiles are shown in figures 4 and 5, respectively. 
In every case the alternate theory predicts more cross 
modulation than the original theory. [When h = 
45 km and N 4 is used in condition 1 (see fig. 4), then 
T & n/T OTig =10; when h o = 30 km and N 5 is used in 
condition 2 (see fig. 5), then T &lt /T orig =5.5.] This 
indicates that less ionization is required in the lower 
regions to produce a given amount of cross modula- 
tion when the alternate theory is used in place of the 
original theory. 

The " crossover" point, mentioned in point 5 of 
the last section provides, at present, the most re- 
liable information that can be obtained from an 
experiment in ionospheric cross modulation. This 
information concerns the electron collision frequency, 
and the next paper is completely devoted to this 
subject [Benson, 1964b]. The total cross modula- 
tion T, as given by the integral in (31), reverses 
sign because the term (dk w /dv m ) in the integrand 
reverses sign. This term, which is given by (7), 



changes sign at the point where a= 



_27rf w ±co L _ 



= 2.18. 



Since j w and o) L =\s cos <£| are known, the electron 
collision frequency v m can be determined at the 
height where (dk w /dv m ) = 0. The location of this 
height can be estimated by observing the value of 
h that causes the complete integral in (31), i.e., the 
observed cross modulation, to change sign. As 
mentioned in the last section, this "crossover" 
point is only slightly dependent on the form of the 
electron density profile when the original theory of 
cross modulation is used. Unfortunately, the alter- 
nate theory of cross modulation is valid only when 
f d <C<Cv m and thus cannot be used in the vicinity 
of, or above, the crossover region. The crossover 
point should not be altered by more than a few km, 
however, even if a more general expression for 
Ad/6, which holds for the entire D region, is used 
since this point is determined mainly by the term 
(dfcu/dO in (31). In support of this last state- 
ment are the observations of this cross-modulation 
"crossover" point that yield values for the electron 
collision frequency that are in agreement with 
other independent observations [Benson, 1964b]. 

In summary, it appears that the general shape of 
the cross-modulation curve can be predicted from 
the original theory of cross modulation [Bailey and 
Martyn, 1934], but that the absolute magnitude of 
the cross modulation is larger than is predicted by 
this theory, and that it is in better agreement with 
the alternate theory proposed by Huxley [1953] in 
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the lower D region. It is apparent that a different 
expression, that will hold for the entire D region, is 
required for the term Ad/6 that enters in (31). This 
implies that a different expression is required for 
the term v m bQ in (22). A recent theoretical in- 
vestigation by Altshuler [19(33] shows that (he 
electron cooling law is determined by inelastic 
electron collisions and that the term v m hQ is pro- 
portional to AQ/Q l/2 rather than to AQ as predicted 
by the original theory of cross modulation or to 
v m (AQ) 2 /Q 2 as predicted by the alternate theory of 
cross modulation. The energy dependence of this 
term in Altshuler's analysis is very similar to the 
energy dependence of this term in the original 
theory; thus, a cross-modulation profile based on 
Altshuler's results should have approximately the 
same general shape as a profile based on the original 
theory of cross modulation. The absolute magni- 
tude of the cross modulation, however, will depend 
on the constant terms (which were not evaluated in 
Altshuler's analysis) in the expression for v m bQ. 
This problem will be considered further in future 
publications concerning the cross-modulation ex- 
periment at College, Alaska. 



I am grateful to W. L. Flock for encouraging this 
work and for many helpful discussions, and to 
G. C. Rumi for his critical correspondence. 

This work was supported by the Air Force Cam- 
bridge Research Laboratories Office of Aerospace 
Research under Contract No. AF 19(604)-8833 and, 
in part, by the Defense Atomic Support Agency, 
Washington, D.C. under Web No. 04.0063. 

7. Appendix 

Huxley's [1953, section 3] alternate theory of cross 
modulation is based on laboratory measurements 
[Crompton, Huxley, and Sutton, 1953] of the temper- 
atures, energy losses, and collisional frequencies of 
electrons drifting through air due to a uniform and 
constant electric field. After converting from cgs to 
mks units and from v to v m (y= 1.128 v m ), to be con- 
sistent with the present paper, his basic equations 
become 



v m =ak T p(mmHg) 

where k T =Q/Q and a=8.30X10 7 , 

bZ(v/m) 



k T =l + 



p(mmHg) 



a) 



(2) 



where Z is the steady applied electric field and 
6=0.33, and 



W m =(lA9XU) n ) ^^m/sec 

Vm 

where W m is the most probable drift velocity. 



(3) 



The average power w supplied to an electron by the 
constant electric field Z is 

w=Z(v/m) e (coulomb) W m (m/sec) j/sec. (4) 

In equilibrium conditions, w= v m dQ, from (22), and the 
above equations can be combined to give 

V m 0y=±? 7^ V m 



Q 2 



where 



^ = (1.49X10")[.(coulomb)1 =318><1 

a z b z 

as stated in (27). 

Next, w (the power supplied to each electron) must 
be expressed in terms of P d (the peak power radiated 
by the antenna of the disturbing transmitter). If 
the antenna has a gain g, then the disturbing power 
per unit area at the height hf is given by 



where 



P<gF(h') 

F{h>) = e-J> idh 



takes into account the absorption of power up to the 
height W . The power per unit area absorbed in an 
infinitesimal layer of thickness clh' is given by 

Thus, the power supplied to each electron »i* - is 

P d (fK d - f h, 2K d dh 

w= — — — P J° 

2irN{h') 2 e 

In (30) the absolute value of w is given since the 
quantity (w) 1/2 enters in the equations for AQ in the 
alternate theorv of cross modulation. 
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